Ultrafast interatomic electronic decay in multiply excited clusters 
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An ultrafast mechanism belonging to the family of interatomic Coulombic decay (ICD) phenomena 
is proposed. When two excited species are present, an ultrafast energy transfer can take place 
bringing one of them to its ground state and ionizing the other one. It is shown that if large 
homoatomic clusters are exposed to an ultrashort and intense laser pulse whose photon energy is in 
resonance with an excitation transition of the cluster constituents, the large majority of ions will be 
produced by this ICD mechanism rather than by two-photon ionization. A related coUective-ICD 
process that is operative in heteroatomic systems is also discussed. 

PACS numbers: 31.70.Hq, 32.80.Rm, 36.40.-c, 32.80.Wr 
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The rapid development during the last decades of very 
intense light sources with extreme short pulse duration 
opened a new era in the study of radiation-matter in- 
teraction. Studying the interaction of intense fields with 
matter brought to the discovery of a whole plethora of 
new physical phenomena, like high-harmonic generation, 
above-threshold ionization, or tunneling ionization, to 
name only a few. In the same time, the progress in gen- 
erating extremely short pulses gave the scientific commu- 
nity a powerful tool to monitor and control the electron 
dynamics in atomic and molecular systems and to study 
processes that take place on a time scale in which the elec- 
tronic motion is still disentangled from the slower nuclear 
dynamics (for recent reviews see, e.g., Refs. [H, [1]). A 
number of free-electron lasers are in operation today pro- 
viding extremely bright, coherent, and ultrashort pulses 
in the VUV regime. Exposed to such highly intense 
pulses, atomic and molecular systems will absorb a large 
amount of photons triggering various dynamical effects. 
In this letter we will restrict ourselves to situations where 
the single-photon energy in the pulse is not high enough 
to directly ionize the system. It is well known that even 
in this case the system can be ionized by a multipho- 
ton ionization mechanism. The multiphoton ionization 
(MPI) results from the ability of quantum systems to 
absorb several and even many photons, whose individ- 
ual energies are insufficient to ionize the system. The 
combined energy of the absorbed photons, though, suf- 
fices to eventually eject one or many electrons from the 
system. During the last decade the MPI has been in- 
tensively studied also in composite systems, like clusters, 
employing the new powerful laser sources (for a review 
see, e.g. Ref. Q). However, little attention was paid to 
other mechanisms that can lead to a multiple ionization 
in an atomic or molecular cluster irradiated by an intense 
laser pulse. 

In this letter we aim at discussing a hitherto unrec- 
ognized mechanism for producing ionized species in ho- 
moatomic or homomolecular clusters exposed to an in- 
tense laser pulse, which in many cases can be by far 
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FIG. 1: Schematic representation of the process, a) Two sub- 
units of the system are excited by absorbing two photons, b) 
One of the constituents of the system de-excites transferring 
the energy to the neighbor which uses it to emit its excited 
electron. 



the dominating one. For simplicity we will consider an 
atomic cluster, but we stress that the process is gen- 
eral and not restricted to atomic systems. Let us take 
a homoatomic cluster and irradiate it with a short and 
intense laser pulse with photon energy below the ioniza- 
tion threshold of the cluster constituents but in resonance 
with one of their excited states. A fraction of the clus- 
ter constituents will be ionized by MPI but, since we are 
at resonance, the large majority of atoms will be exited. 
The well known effect of Coulomb blockade will not play 
a significant role here, since we suppose that the system 
is exposed to a very short pulse, i.e. to a broadband ex- 
citation. Thus, there will be many excited atoms in the 
cluster whose neighbor is also excited. Having two ex- 
cited atoms in close proxmity the following interatomic 
electronic decay mechanism is conceivable. One of the 
atoms is de-excited, the energy is transferred to the other 
one which uses it to emit its excited electron. Thus, at 
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the end of the process, one of the atoms is ionized and the 
other one has returned to its ground state. The process 
is pictorially represented in Fig. [T] and can be written in 
short as: 

A* ---A* ^ A---A+ + e-. 

This process bears similarities with the interatomic (in- 
termolecular) Coulombic decay (ICD) predicted theoret- 
ically more than ten years ago [J| and since then stud- 
ied very intensively both theoretically and experimen- 
tally (see, e.g. Refs. The ICD is a very efficient 
electronic decay mode of inner-valence ionized atoms or 
molecules embedded in an environment. Inner- valence 
ionized states usually have energies below the double ion- 
ization threshold and, thus, cannot autoionize. However, 
here the environment plays a critical role. When the 
initially ionized atom or molecule has neighbors, like in 
a cluster, an electron from a higher level may fill the 
vacancy and the released energy can be transferred to 
a neighbor form which a secondary electron is emitted. 
Thus, the creation of a single hole in one of the subunits 
in the system leads to the formation of two positively 
charged subunits that repel each other typically leading 
to a Coulomb explosion that disintegrates the system. 
The process is ultrafast with typical lifetimes of few to 
few tens of femtoseconds, quenching all other energeti- 
cally allowed relaxation modes of the system. The dis- 
covery of the ICD revealed a whole zoo of related phe- 
nomena, involving both energy and electron transfer and 
initiated by single or multiple ionization, as well as by 
inner- or outer valence excitation (for recent review, see 
Ref. 0). Although these processes have different names 
and acronyms, we will refer here to all these phenomena 
as ICD in order to make the text more transparent. The 
only, but important difference of the process proposed 
here and the ICD phenomena studied until now is that 
the ICD assumes an excited system interacting with a 
non-excited environment, while in the process sketched 
in Fig.[T]the distinction between system and environment 
is not possible. On the contrary, both constituents are 
equally suitable to undergo an electronic decay. However, 
we will refrain from giving a new name to the process dis- 
cussed in this letter and will refer to it as ICD. 

The important question is, of course, whether this ICD 
process is efficient enough and can compete with the 
other possible de-excitation modes (e.g. photon emis- 
sion) in the dimer or, even more interesting, in a large 
cluster. To estimate that we have to calculate the rate 
of the process, or the decay width F. The easiest way 
to estimate the decay width is to consider the process 
within the simplified but insightful picture of interaction 
between two dipoles via a virtual-photon exchange. A 
virtual photon is emitted as a result of the de-excitation 
of one of the excited atoms and then absorbed by the 
other excited atom causing its ionization. The virtual 



photon exchange picture, which is correct at large in- 
teratomic distances, enables the derivation of analytical 
formulae for the decay width ll| . Such formulae exhibit 
dependence {R being the interatomic separation) 
with a prefactor specific for the emitting and absorbing 
constituents and accounting for the dipole selection rules 
of the involved transitions. The derivation of such an ex- 
pression is straightforward using the procedure explained 
in detail in Ref. [l^l . 

In order to illustrate the efficiency of the ICD mecha- 
nism proposed here we consider a concrete example. Let 
us take a neon dimer in which both of the neon atoms 
are m their first excited state Ne*(2p-i3s)-Ne*(2p-i3s). 
The energy of Ne*(2p-i3s) is about 16.7 eV above the 
ground state while the ionization potential (IP) of the 
neon atom is about 21.6 eV. Thus, the 3s — >■ 2p tran- 
sition in one of the neons will release enough energy to 
ionize the other one, emitting an electron with kinetic en- 
ergy of about 11.8 eV. Averaging over the multiplicities of 
the initial states and summing over the final states we ob- 
tain for the total decay width of the system Ne*(2p~^3s)- 
Ne*(2p~^3s) as a function of the internuclear distance R 
the following expression (in atomic units) 



T{R) = 



3c fa 1 



(1) 



where / is the oscillator strength of the 3s ^ 2p transi- 
tion, cr is the ionization cross section of Ne*(2p~^3s), c 
is the speed of light, and w is the virtual-photon energy. 

The values of the quantities entering Eq. ([T]) are known 
from the literature - the oscillator strength for the 3s — > 
2p transition in neon is 0.16 [l3l | and the photoionization 
cross section of Ne*(2p-i3s) is about 0.18 Mb Q. At 
the Ne2 equilibrium distance of 3.1 A the decay width 
for the process is 0.24 meV, which implies a life time of 
about 2.8 ps. This is 3 orders of magnitude faster than 
photon emission, which is known to be about 2 ns [isj . 
and thus ICD is by far the dominant relaxation pathway 
in the dimer. The results obtaind by the virtual-photon 
model are correct at large interatomic distances R. They 
are very promissing, in particular, since it is known from 
previous studies [llj that such kind of asymptotic for- 
mulae underestimate the decay rates around equilibrium 
distances due to neglecting the orbital overlap. 

In order to have more reliable values for the decay rate, 
we used the ah initio method, known as Fano-Stieltjes 
approach [l^. In this method the boundlike and the 
continuumlike components of the wave function of the 
decaying state are constructed using the Green's func- 
tion formalism, and the problem of the normalization of 
the continuum wave function is addressed by using the 
Sticltjes imaging technique (see Ref. [l^ for details). 
The ab initio results are shown in Fig. [2] together with 
the predictions of the virtual- photon model, Eq. ([1]). For 
a reference, the atomic fiuorescence decay width is also 
shown in the figure. We see that up to about 9 A of in- 
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FIG. 2: Total ICD width T for the system Ne* {2p-^3s)- 
Ne*(2p-^3s) compared to the prediction of virtual photon 
model, Eq. IT}. The atomic fluorescence decay width is indi- 
cated by a horisontal line, while the equilibrium interatomic 
separation by a vertical one. Note the double logarithmic 
scale used. 



ternuclear separation, i.e. about 3 times the equilibrium 
distance, the asymptotic formula largely underestimates 
the ICD decay width. When the width becomes very 
small (i.e., at large internuclear separation) the ab initio 
method suffers fron numerical instabilities and cannot be 
safely employed. It is also around 9 A distance where the 
radiative decay becomes competitive. At the equilibrium 
distance of the neon dimer the ab initio computation pre- 
dicts a total decay width of 5.4 meV which is more than 
20 times larger than the virtual photon result. This decay 
width corresponds to a life time as short as 122 fs, which 
means that the ICD sets in before the nuclear dynamics 
play a role. 

Let us now comment on larger clusters. Most impor- 
tantly, since the total decay width is a sum of the partial 
widths of all possible decay channels, it is clear that if 
we have more than two interacting excited atoms the 
ICD process will become dramatically faster [HI, [l3| . In 
(Ne*)4, for example, there are 12 open channels, which 
suggests that the ICD life time in this cluster will be 6 
times shorter than that for (Ne*)2. Thus, in big clus- 
ters, where a resonant intense laser pulse will produce a 
large number of excited atoms, the ICD mechanism will 
be extremely efficient. 

Once we have seen that the ICD process is ultrafast 
and can be expected to outperform other possible ways 
of relaxation, let us return to the question of the compe- 
tition between the ICD and the MPI in the production 
of positive ions in a cluster irradiated by a laser pulse 
with high density of photons. For that purpose, it is il- 
luminating to consider again a concrete example. Let a 
Neiooo cluster be exposed to a short and intense laser 
pulse with photon energy of 16.7 eV, i.e. resonant to the 



2p — > 3s excitation of the neon atom. We can estimate 
the number of excited atoms and the number of those 
ionized by two-photon ionization in the cluster after the 
pulse by solving the following system of rate equations 



dN{t) 
dt 



-aoHt)N{t)-a2^^it)N{t), 



= <yomN{t) - a,<f{t)N^*\t), 

dN(+\t) 



dt 



CTl 



$(t)iV(*)(i)+(T2$2(t)iV(t). (2) 



In Eqs. © N{t), N^*\t), and N^+\t) are the number 
of neutral, excited, and ionized by two-photon ionization 
atoms as a function of time, respectively, while cto is the 
absorption cross section, ai is the photoionization cross 
section of Ne*{2p~^3s), a2 is the two-photon ionization 
cross section, and denotes the photon flux which 
contains the information on the temporal profile of the 
pulse. In order to obtain quantitative results, one has to 
consider also the spatial profile of the pulse and the ge- 
ometry of the irradiated cluster. However, we aim here 
at making only an estimate of the ratio between A''^*) 
and N'^'^^ after the pulse and that is why we will use 
the rather simplified picture of a rectangular pulse with 
intensity 10^^ W/cm^ and duration 50 fs, ignoring the 
dependence of the laser-cluster interaction on the spa- 
tial profile of the pulse and the geometry of the cluster. 
In this case, Eqs. ^ can be easily solved and using the 
atomic data, ctq ~ 273 Mb [H, ai « 0.18 Mb [li|, and 
(72 ~ 2 X 10~^^ cm^'s [l3|, one obtains that in the Neiooo 
cluster after the pulse 991 atoms will be excited and only 
3 will be ionized by a two-photon ionization. Since, as we 
saw, the ICD process is very efficient, one would expect 
that every pair of Ne*(2p~^3s) will undergo ICD produc- 
ing about 495 neon ions. Thus, the ratio of the neon ions 
produced by ICD and those produced by a two-photon 
ionization is about 166:1. It is clear that by increasing 
the laser intensity one will produce more ions by two- 
photon ionization, while decreasing it will favor the ion 
production via ICD mechanism. For example, with a 
laser intensity of 10^^'^ W/cm^ this ratio is 13.7:1. We 
see that even at these relatively high intensities, the ICD 
mechanism is still by far the dominant source of ionized 
species in the cluster. It is clear that even at higher inten- 
sities, the production of ions via ICD has to be taken into 
account when interpreting experimental results. Indeed, 
the peak intensity is achieved only in the focal point of 
the laser which usually is much smaller than the interac- 
tion region. A large fraction of the clusters, thereby, will 
be exposed to a less intense field where the ICD is the 
dominant ion-production source. 

At the end we would like to comment briefly on an- 
other possibility to create ionized species in multiply ex- 
cited clusters which will be operative in the case of het- 
eroatomic systems. In the case when the de-excitation 
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FIG. 3: Schematic representation of the collective process dis- 
cussed in the text for heteroatomic clusters. Two constituents 
of the system de-excite simultaneously transferring the energy 
to a third species and ionizing it. 



Moreover, we showed that if large clusters are exposed 
to an ultrashort and intensive laser pulse (10^^ — 10^'^ 
W/cm^ in the present example) which is in resonance 
with an excitation transition of the cluster constituents, 
the large majority of ions will be produced by this ICD 
mechanism rather than by two-photon ionization. In ad- 
dition, we proposed a collective ICD process that can 
take place in heteroatomic or heteromolecular systems 
also yielding ionized species. We hope that our work will 
trigger more theoretical and experimental investigations 
of these ICD effects in systems exposed to ultrafast laser 
pulses with high density of photons. 
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(FP7/2007-2013) / ERC Advanced Investigator Grant 
n° 227597. 



energy of an excited atom is insufhcient to ionize an- 
other atom, a process related to the recently discussed 
coUective-ICD [19j can take place. In the collective-ICD 
process two inner- valence ionized species de-excite simul- 
taneously transferring their "collective" energy to a third 
neighbor and ionizing it. In analogy, one can think about 
a collective-ICD where two excited atoms or molecules 
de-excite simultaneously and the released energy is used 
by a third atom or molecule to eject one of its electrons, 
see Fig. |31 An important point to note is that, in contrast 
to the former case, in the case of collective-ICD from ex- 
cited species the process will not have to compete with 
the Coulomb explosion dynamics of the two neighbor- 
ing ions. It is clear that the collective-ICD from excited 
species will be energetically open when 2E{A*) > IP{B). 
This implies that A should be different from B since we 
supposed that the ICD process of Fig. [T]is energetically 
closed. Of course, if the ICD channel is open, the collec- 
tive decay can also take place, but since it involves three 
electrons, its importance compared to the two-electron 
ICD process will be low. 

Let us conclude. In this letter we proposed a hitherto 
unrecognized mechanism for producing ionized species 
in multiply excited atomic or molecular clusters. The 
mechanism belongs to the family of interatomic (inter- 
molecular) Coulombic decay phenomena and consists of 
an ultrafast energy transfer between two excited species, 
bringing one of them to its ground state and ionizing 
the other. We showed that the process is ultrafast (in 
the femtosecond time regime) and as such is extremely 
efficient compared to other possible relaxation modes. 
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